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Abstract. I examine the uncertainty of perturbative QCD factorization for (un)polarized hadron 
structure functions in deep inelastic scattering at a large value of the Bjorken variable xb- The 
focus will be on Target Mass Corrections and Jet Mass Corrections in the collinear factorization 
framework. 
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Precise partem distribution functions (PDFs) at large parton fractional momentum x 
are vital for understanding the nonperturbative structure of the nucleon, and the effects 
of color confinement on its partonic constituents. For instance, the d/u quark distri- 
bution ratio near x = 1 is very sensitive to the nature of the quark-quark forces in the 
nucleon, as are parton-parton correlations in the nucleon wave function. Moreover, the 
magnitude and sign of the ratios of spin-polarized to spin-averaged PDFs, Au/u and par- 
ticularly Ad I d, in the limit x — > 1 reflect the nonperturbative quark-gluon dynamics in 
the nucleon, and can shed light on the origin of the nucleon's spin. Improved PDFs at 
large x will be crucial at the upcoming 12 GeV Jefferson Lab Upgrade for the planning 
and analysis of a host of experiments, as well as at upcoming neutrino oscillation ex- 
periments such as MINOS and MINERvA at Fermilab to control and reduce systematic 
errors to a minimum. They will also directly impact experiments at hadron colliders such 
as the Large Hadron Collider at CERN, by facilitating accurate extraction of new physics 
signals as excess compared to benchmark perturbative QCD (pQCD) calculations, and 
by better constraining PDFs at lower x and larger Q 2 via pQCD evolution. 

Parton distribution functions can be extracted from experimental data through global 
QCD fits, which combine data from many different processes and observables, and an- 
alyze them by means of pQCD calculations. Currently, the unpolarized parton distri- 
butions are well-determined for x < 0.5 for quarks and x < 0.3 for gluons. To better 
constrain them at large x it is necessary to study hard scattering processes near kine- 
matic thresholds, such as Deep Inelastic Scattering (DIS) at large Bjorken invariant xb 
and low 4-momentum transfer squared Q 2 , or Drell-Yan (DY) lepton pair production 
at large Feynman variable xf. The large-* and low-<2 2 region is even more important in 
the case of parton helicity and transversity distributions, which measure the parton's spin 
distribution and correlation to its orbital angular momentum, since much of the presently 
available experimental data lie in this kinematic domain. 

Theoretical improvements in calculating leading-twist perturbative cross sections are 
needed in order to confidently access large-* / low-<2 2 PDFs. In this talk, I use the 
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collinear factorization framework to formulate and quantify corrections to DIS structure 
functions arising from target and jet mass corrections |Q1 0]. Threshold resummation, 
higher twists, quark-hadron duality in the pre-asymptotic domain, and nuclear effects 
such as binding, Fermi motion and nucleon off-shellness are left for future efforts. 

Target mass corrections (TMC) for spin-averaged nucleon structure functions were 
first considered by Georgi & Politzer within the Operator Product Expansion (OPE) [0, 
0]. For spin-dependent scattering they were evaluated within the same OPE formalism 
in Refs. Q4J, |5Q. One of the limitations of the OPE formulation of TMCs is the so-called 
"threshold problem", in which the target mass corrected structure functions remain 
nonzero at xb > 1- This arises from the failure to consistently incorporate the elastic 
threshold in moments of structure functions at finite Q 2 , resulting in nonuniformity of 
the Q 2 — > °° and n — > °° limits, where n is the rank of the moment. As a consequence, 
after performing an inverse Mellin transform on the moments, the structure functions 
acquire incorrect support at large xb la]. Various prescriptions have been considered to 
tame the unphysical behavior as xb — > 1 • However, these approaches are not unique and 
sometimes introduce additional complications. 

An alternative approach, which avoids the threshold ambiguities from the outset, in- 
volves formulating TMCs directly in momentum space using the collinear factorization 
(CF) formalism yj]. A further advantage of this formalism is that it can be readily ex- 
tended to processes such as semi-inclusive DIS, where the OPE is not available, and 
indeed to any other hard scattering process. Additional corrections to structure functions 
at large xb, such as discussed in the introduction, can also be naturally incorporated 
together with TMCs. At any order in the strong coupling constant a s , and considering 
only massless partons, the CF target-mass corrected transverse and longitudinal structure 
functions, Ft and Ft, read 

rU x B rfr f£ \ 
F t ,l(xb,Q 2 ,M 2 )=£ ^ ™h flT>L i^-,Q 2 J cp f (x,Q 2 ) , (1) 

where / is the parton flavor, M the target's mass, £ the Nachtmann scaling variable, 

K = - 2x ^- 2 r 2 = 44m 2 /2 2 , (2) 

hf the perturbatively calculable coefficients, and <j)y the unpolarized PDFs. The polarized 
gi structure function reads 

gl {x B ,Q\M 2 ) = T — 2 £ jf -SU (^Q 2 ) A<P/(^2 2 ) , O) 

where g\j are the perturbatively calculable coefficients, and A<py are the parton he- 
licity distributions. The upper limit of integration x < E,/xb is imposed by combining 
four-momentum and baryon number conservation in the kinematics of the handbag dia- 
gram. As a result, the structure functions are zero when xb > 1, which solves the thresh- 
old problem. Note also that contrary to previous claims 10], Ft,l(xb) ^ FjI(%), and 

S\{ x b) 7^ (1 + y 2 ) (£)» where the superscript (0) indicates "massless" structure 
functions computed with M = 0, and the upper limit on the x integration is x < 1 . 




FIGURE 1. Square plots: Comparison of prescriptions for NLO target mass corrections to Fj_ and 
R = Fi/F\, computed with MRST2002 parton distributions 1 1]. (The "GP" curves are computed in OPE; 
the "nv" curves are computed using 1 as upper bound for the x integration in Eq. ([T). Vertical plots: Same 
for the LO g\ structure function (left) and Ai polarization asymmetry (right) 121. For Ai, the shaded band 
for the OPE result indicates the uncertainty due to the use of the Wandzura-Wilczek relation (lower bound) 
or the identity g \ + g2 — (upper bound) in deriving A i . 



An estimate of TMC in collinear factorization for polarized and unpolarized struc- 
ture functions, and for the A\ = (1 + y 2 )g\/F\ asymmetry is presented in Fig. [H and 
compared to the TMC computed in OPE. The numerical results for the target mass cor- 
rections to the F2 and g\ structure functions are qualitatively similar in the CF and OPE 
computations. At leading order (LO), the differences of up to 30% for low Q 2 values dis- 
appear for Q 2 > 10 GeV 2 . At next-to-leading order (NLO), the differences are generally 
larger than at LO for intermediate xg values. In all cases the TMCs remain significant 
at %b > 0.7 even for Q 2 > 10 GeV 2 , and need to be taken into account when analyzing 
large-jcs data. In the longitudinal Fi structure function the CF and OPE corrections are 
larger by themselves, and differ by more than factors of several at large and not-so-large 
xb, even at large Q 2 values usually considered TMC-free. The choice of TMC scheme 
will affect, e.g., the determination of the gluon PDF, to which Fi is much more sensitive 
than F2, and is not a negligible matter. Since the TMCs are qualitatively similar for the 
gi and F\ structure functions, they largely cancel in the A\ asymmetry, although the sign 
of the correction is opposite in the CF and OPE approaches over most of the range of xb, 
but can be as large as 20% at low Q 2 , disappearing at Q 2 > 10 GeV 2 . The commonly 
used approximation relating A 1 directly to the longitudinal asymmetry Ay can be shown 



to be equivalent to assuming A \ « A\ yD. Fig. Q] shows that it breaks down at low Q 2 : 
an accurate measurement of g\, e.g., at JLab, requires both Ay and A± asymmetries. 

Even after solving the threshold problem as described, the structure functions remain 
unphysically positive as xg^l because the assumption that at LO the final state current 
jet is made of a single massless quark makes the LO perturbative coefficients propor- 
tional to 8(x — <fj /xb). Such an assumption is clearly unphysical because the quark has 
to hadronize due to color confinement, so that the current jet will have an invariant mass 
m 2 . Using this mass in the kinematics of the LO handbag diagram, and assuming that 
the invariant jet mass has a probability distribution J m , one obtains 

F{ MC (x B ,Q\M 2 ) = j^ Q2 dm)j m {m))F^^{\ +m 2 J /Q 2 ),Q 2 ) . (4) 

If the "jet function" J m (m 2 ) is a sufficiently smooth function of m 2 , one obtains 
Fj. mc (xb, Q 2 ,M 2 )^0 as xg— >1, as it should be. The jet mass corrections (JMC) so in- 
troduced are of order 0(m 2 /Q 2 ), and it is easy to see that in the limit Q 2 ^> (m 2 ) one 
recovers the structure functions with TMC only. It turns out that the ansatz © is theo- 
retically correct, and that the jet function is the spectral function, J2, of a vacuum quark 
propagator, 

Jdm 2 J 2 (m 2 )27td(l 2 -mj)e(l ) = J d 4 ze lzl Tr[r (0|^(z)y(0)|0)] , (5) 

additionally smeared by soft momentum exchanges with the target jet [1]. A toy model 
estimate of these Jet Mass Corrections shows that they may in fact be comparable in size 
to the TMC. However their phenomenology is just at its beginning. 
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